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(— i . Abstract. The steady states of two gases of hard spheres or disks separated by an 

O | adiabatic piston in presence of a temperature gradient are discussed. The temperature 

field is generated by two thermal walls at different temperatures, each of them in 
contact with one of the gases. The presence of the piston strongly affects the 
hydrodynamic fields, inducing a jump in its vicinity. A simple kinetic theory model 
is formulated. Its predictions are shown to be in good agreement with molecular 
dynamics simulation results. The applicability of the minimum entropy production 
principle is analyzed, and it is found that it only provides an accurate description of 
the system in the limit of a small temperature gradient. 
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The adiabatic piston consists of a container filled with a gas that is divided into 
two compartments by a freely moving adiabatic piston pp. Usually, the system is 
isolated and initially prepared with the two gases in both compartments at different 
independent states. Then, the relaxation to equilibrium is studied. It is observed 
that the system first converges towards a state of mechanical equilibrium with both 
gases having the same pressure. Then, the pressure fluctuations that are asymmetrical 
because the temperatures of the gases are different, drive the system very slowly to 
thermal equilibrium (2j |3]. Of course, the situation is much more complex and less 
understood when the system can not relax to equilibrium due to some imposed external 
conditions. Very recently, several configurations consisting of granular gases separated 
by an adiabatic piston have been considered [H [5j [6] . Granular gases are inherent non- 
equilibrium systems, since no equilibrium state is possible for them and, for this reason, 
they have been intensively used to investigate many fundamental issues in the context 
of non-equilibrium statistical mechanics [7J|S]. Nevertheless, it must be realized that the 
irreversibility in granular gases has an internal origin, and it is not related, in principle, 
to some imposed boundary or initial conditions, as it is the case in molecular systems. 
Actually, as a consequence of the above, gradients and inelasticity are related in steady 
granular gases and this leads to a quite peculiar rheological behavior. 

The aim of this paper is to study the steady state of two molecular gases separated 
by an adiabatic piston, when each of the gases is in contact with a thermal wall at 
different temperature. Attention will be restricted to the final state of the system, 
whithout studying the way in which it is reached. The questions addressed are related 
with the modification of the hydrodynamic steady profiles because of the presence of 
the piston and the steady position of the latter. 

The system considered is a rectangular (d = 2) or cylindrical (d = 3) container of 
length L x . It is divided into two compartments by a movable adiabatic piston of mass 
M, constrained to remain perpendicular to the axis of the system, taken as x-axis. By 
adiabatic it is meant, as usual, that the piston has no internal degrees of freedom and, 
therefore, it can not transmit energy when it is at rest. A sketch of the system is given 
in Fig. [TJ In each of the two compartment there are N elastic hard disks (d = 2) or 
spheres (d = 3) of mass m and diameter a. Collisions between particles and the piston 
as well as collisions of the particles with the lateral walls of the container are also elastic. 
Moreover, the motion of the piston occurs without friction with the lateral walls. 

The walls of the container located at x = and x = L x and perpendicular to the x- 
axis are thermal walls with temperatures T 10 and T 2 o, respectively. At the microscopic 
level of description, a thermal wall is modeled by assuming that whenever a particle 
collides with it, a new velocity is assigned to the particle drawn from a Maxwellian 
distribution with a second moment defined by the temperature of the wall [9j [TO] . 
More specifically, the velocity distribution of the particles leaving a thermal wall with 
temperature T i0 , i — 1, 2, is given by 

PM = e [(-ir V] (^f *«p (-^) , (i) 
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Figure 1. Sketch of the system considered in this work. The left and right walls are 
thermal with temperatures T\q and T20, respectively. 



where 6(x) is the Heaviside step function and ks is the Boltzmann constant. At the 
macroscopic level, it will be assumed that the fluid in the vicinity of a thermal wall has 
the same temperature as the wall. 

In the following, attention will be restricted to steady states in which the gases 
at both sides of the piston are very dilute, so that they verify the local equation of 
state pi = UikBTi, where pi, Hi, and T i} are the local pressure, number density, and 
temperature of the gas in compartment i, i = 1, 2. Also, in the steady states considered 
only gradients of the hydrodynamic fields in the x direction are present. 

The hydrodynamic Navier-Stokes equation for the velocity field, when applied to 
the above steady states, implies that the pressure must be uniform in each container. In 
addition, mechanical equilibrium of the piston requires that the pressure on both faces 
be the same, 

Pl=P2= P- (2) 
The energy balance in each container reads 

where q X i is the heat flux along the x direction in compartment i. According to Fourier 
law, 

dT- 

q xi = -K(Ti)-^ . (4) 
The thermal conductivity k(T) in the dilute limit being considered is 

From Eqs. ©-© it follows that 

f)2rn3/2 

-a£- = o, (6) 
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so that 

Xf /2 (x) = aiX + k, (7) 

with cij, bi constants. But stationarity also requires that the heat flux be the same in 
both compartments, since otherwise the energy of the piston could not be constant. 
Then, q x i = q x2 and using the Fourier law this leads to 

ai = a 2 = a. (8) 

Now the boundary conditions at the thermal walls, Xi(0) = Tio and T2(L X ) = T 2 q, are 
imposed, to get 

T^(x) = ax + T^ 2 , (9) 
T 2 3/2 (x) = -a(L x .-x)+T 2 3 /2 . (10) 



Equation ([9]) holds for < x < xp, and Eq. ( UOl for xp < x < L x , where xp is the 
steady position of the piston. 

The boundary conditions associated to the piston must also be introduced. Then, 
expressions for the energy flux between the piston and the two gases are needed. This 
flux appears as a consequence of the velocity fluctuations of the piston [TTj . Assuming 
that they are Gaussian in the steady state, something that is confirmed with very good 
accuracy by the numerical simulations, in the limit m/M <C 1 it is obtained that the 
energy flux from the piston to the gas in compartment i is given by [5] 

Q^-2 f^) 1/2 T.M-T P ^ 

V t ) Tl /2 {x P ) M' 1 1 

where Tp is the temperature parameter of the piston, defined from the second moment 
of its velocity distribution, and the temperature Ti(xp) is to be understood as the 
temperature of the gas in compartment i in the limit x — > xp taken from inside the gas. 
Stationarity of the piston yields Qi + Q2 = 0, i.e. 

Ti(x P ) - T P T 2 {x P ) - T P 



Tl l2 {xp) T 2 1/2 (x P ) 
or 

T P = [T^xp^xp)] 1 ' 2 . (13) 

A second condition at the piston follows from the continuity of the energy flux, implying 
that 

q x i{x P ) = -Qi- (14) 
Use of Eqs. ® and ^ yields 

T l 2 /2 {x P )-Tl l2 {xp) = C-, (15) 

p 



(12) 



where 



d-2 



c _ d(d + 2) 2 T(d/2)n- — k B M 
48^2(d - l)ma d - 1 
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The last needed condition is that the number of particles in each compartment is fixed 
to N, so it must be 

Sp / dxni(x) = Sp dxn 2 (x) = N, (17) 

JO Jxp 

where Sp is the area (d = 3) or length (d = 2) of the piston and, therefore, also the 
section of the container. By employing the local equation of state and Eqs. ([9]) and (fTUj) . 
the above conditions (1T51) and ( fTTl) are seen to be equivalent to 



T 1 1/2 (xp)+T 2 1/2 (xp)=T 1 1 /2 + T 2 1 /2 



and 



with 



Tl /2 {xp) = T l l ' 2 + D- ) (19) 
V 

Nk B , . 

B = W < 20 > 

Equations (EES]), ([18]) . and ffl9]) form a closed set of equations for the unknown Ti(xp), 
T 2 (xp), and a/p. Solving it gives 

rpl/2 rpl/2 

Tl'\x P ) = Ttf + D T ™ D+ T g , (21) 

rpl/2 _ rpl/2 

T} / \xp) = Ttf-D 1 ™ D+ 1 g , (22) 



1/2 _ ^1/2 

(23) 



20 — J 10 



p 2D + C 

The only remaining task is to identify a (or p) and the steady position of the piston x P . 
By means of Eqs. ([§]) and fllOp it is found 

T^ 2 (xp)~T^(xp)-T^ 2 + T^^ 



x P 



Tt ,2 (x P )-T- 



3/2 
10 



Tt /2 (x P ) - T* /2 (x P ) - Tiq 2 + T 2 3 /2 ' (25) 

Therefore, as a consequence of the presence of the piston the slope a of the temperature 
profile decreases as compared with its values in a gas without piston, being the same 
the temperatures of the two thermal walls. 

To check the accuracy of the above theoretical predictions, molecular dynamics 
(MD) simulations of a system of hard disks (d = 2) have been performed, using an 
event driven algorithm [12]. The number of particles in each of the two compartments 
has been iV = 100 and the size of the system is L x = 2Sp = 100cr. The values of 
M/m and T 2 o/Tio have been varied as indicated below. All the simulations started 
with the piston located in the middle of the system (x P = L x /2) and the particles 
homogeneously distributed in each compartment. After a short transient time, of the 
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Figure 2. Steady temperature (blue crosses and solid line) and density (red circles 
and solid line) profiles. The symbols are MD simulation results for a tow-dimensional 
system of 2N = 200 hard disks with L x = 2S P = 100cr, M/m = 10, and T 20 /T w = 1-5. 
The solid lines are the theoretical predictions derived in the main text. 



order of a few collisions per particle, a steady state was always reached. Then, density 
and temperature profiles of the gases in both compartments as well as the position 
and velocity distribution of the piston were measured. To identify the hydrodynamic 
profiles, the system was divided in 20 layers of the same width, parallel to the piston. 
Moreover, the average quantities have been averaged on time (typically 200 registers) 
and also over a given number of trajectories (typically 50). 

In Fig. El the steady temperature and density profiles are plotted for a system with 
M/m = 10 and T 2 o/T 10 = 1.5. The profiles for a system with the same value of the 
mass ratio M/m, but with a much larger temperature ratio, T 2 o/Ti = 5, are given in 
Fig. [3j The symbols are MD results, while the solid lines are the theoretical predictions 
obtained above. The vertical line indicates the average position of the piston. A fairly 
good agreement is observed. It is important to realize that the hydrodynamic fields in 
the vicinity of the average position of the piston are much influenced by the fluctuations 
of the latter, that are not accounted for in the simple model developed here. 

The results for the dependence of the steady position of the piston xp on the mass 
ratio M/ m are shown in Fig HI Three different values of the ratio of the temperatures 
of both thermal wall have been considered, as indicated in the inset. As the mass 
ratio increases the position of the piston tends to a constant value. A similar behavior 
is observed for the temperature of the piston, shown in Fig. |5j In both cases, there 
is again a fairly good agreement between the simulation results and the theoretical 
prediction from the model proposed here. 
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Figure 3. The same as in Fig. [21 but now T20/T10 = 5 




1 1 1 1 

20 40 60 

M/m 

Figure 4. Steady position of the piston xp as a function of the mass of the piston 
M. Three values of the ratio of the temperatures of the two thermal walls have been 
considered, as indicated. The values of L x and TV are the same as in Fig. [21 The 
symbols are MD simulation results while the solid lines are the predictions derived in 
the main text. 
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Figure 5. The same as in FigUJ but now the plotted quantity is the temperature of 
the piston divided by the temperature of the cold wall. 



It seems worth to investigate whether the steady state reached by the system 
considered here can be described by means of a minimum entropy production rule. 
The entropy production function P of the system is p3J [H] 

dT{\x) 



P = S P dxq xl (x 
Jo 



-Qi 



+ Q 2 



Xp 

1 

Tp~ 



dxq xl (x 
1 



dx 

dT 2 -\x) 
dx 



1 



Ti(xp) 
1 

T~ P 



Mx p) T pV ^ 

where the possible discontinuity of the temperature at the piston has been taken into 
account. In this expression, the relations given in Eqs. ([2]), (fiU]) . and (TT7I) are 
assumed to be known since they correspond to boundary conditions. Moreover, the 
limit of small temperature gradients will be considered in the following. Define 

m ^10 + „ T20 — T10 







IT 



(27) 



and introduce the scaled temperature deviations r 1; r 2 , and rp by 



T^xp) = T m + 6t u T 2 (x P ) = T m + 9t : 



2- 



(2f 



Tp = T m + 9t p . 



(29) 
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When these expressions are substituted into Eq. ( 12 6 p and an expansion in powers of 9 
is carried out, it is obtained that 

p / , 6(2mk B y/ 2 S P 

P{Tl,T 2 ,T 3 ) = —— 372" 

TX X l 2 ML x Tm 

x [2D [r 2 + r 2 2 - 2 (n + r 2 ) t p + 2r| 

+ C [rl + r| + 2T m (n - r 2 + T m )] } # 2 

+ £>(# 4 ). (30) 

The values of n, r 2 and T3 minimizing this expression for small namely neglecting 
terms of order # 4 and higher, are 

Tl = -cq^D' T2 = _Tl ' Tp = - (31) 

Using these values, it is easily seen that 

L x ( C + D \ 2Nk B T m 

ai=a2= C + 2D L x ■ (33) 

These results agree with the lowest order expansion in 9 of Eqs. (!23l -( j25l) . as it can 
be easily verified. On the other hand, it is also easily seen that the condition that 
the entropy production P is minimum is not equivalent to the hydrodynamic theory 
developed here outside of the limit of small temperature gradients, leading to results 
that strongly disagree with the numerical simulations. 

In summary, it has been shown that the presence of the adiabatic piston introduces 
a strong discontinuity of the hydrodynamic profiles, increasing their slope in both 
compartments. Moreover, the simple model presented here based on the energy 
balance, reproduces quite well the molecular dynamics simulation results. A minimum 
entropy production requirement only holds in the limit of small temperature differences. 
Although we are not aware of any realization of a situation similar to the one considered 
here, the effects discussed should be observable in experiments, by extending the devices 
used to investigated the relaxation to equilibrium [15j [HI IT7] . 
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